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THERMAL RADIATION FROM A GAS NOT IN LOCAL 


THERMODYNAMIC EQUILIBRIUM 
by John R. Howell 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

The validity of the assumption that local thermodynamic equilibrium is 
present in a gas is examined with respect to its effect on the calculated 
emission of thermal radiation. The spectral distribution of energy emitted 
by the gas is determined when there is departure from equilibrium. Attention 
is confined to the case of insufficient interparticle collisions being avail- 
able to redistribute the energy of absorbed radiation. Thus the equilibrium 
distribution of energy states in the gas cannot be maintained. Emission of 
radiation from such a gas is not in the usual spectral energy distribution 
determined by the equilibrium spectral absorption coefficient and the Planck 
distribution. 

Specifically, the case of the grey gas is used to derive some approximate 
relations for characterizing the conditions where departure from equilibrium 
might be expected. 

INTRODUCTION 

Solutions to radiative energy transfer problems involving absorbing- 
emitting media almost invariably contain the assumption of "local thermo- 
dynamic equilibrium" (LTE). Basically this means that any small gas volume 
is assumed to be in thermodynamic equilibrium so that its temperature can be 


defined. 
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When this assumption is valid, it greatly simplifies radiant transfer 
problems . The emission of energy from a volume element can be specified in 
terms of the Planck spectral energy distribution and the radiative absorption 
coefficient regardless of the spectrum of the flux passing through the element. 
If LTE cannot be assumed, the Planck distribution is not valid and the actual 
distribution, along with the gas absorption coefficient, becomes a function of 
the population of the energy states in the gas. Detailed quantum-statistical 
calculations become necessary to describe accurately the emission and absorption 
of radiant energy. In the general case involving surface-gas and gas -gas 
radiant energy exchange, consideration of such detailed phenomena tax present 
analytical methods past their capabilities. 

It is possible to envision physical systems where nonequilibrium effects 
could be significant. Some of these have been examined by Kulander [l] and 
Oxenius [2]. For example, if very large temperature gradients are present 
such as those in shocks or those expected in proposed gaseous-core nuclear 
rockets, then species migrating to a given position may come from regions at 
widely differing conditions. One temperature could hardly be expected to 
characterize local conditions in such a case. Ferrari and Clark [3,4] have 
examined this situation with reference to hypersonic shock phenomena. 

A less common example is the presence of transients in the incident ra- 
diation. If the transients are quite rapid so that significant changes in 
incident radiation intensity occur over periods of time of the same order as 
those governing the emission process in the gas (electron transitions, vibra- 
tional transitions, etc.), then the emission of energy may greatly lag the 
absorption process. This is a nonequilibrium situation and might be encountered 
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in connection with thermal emission from a nuclear fireball. 

Another case is that of a rarefied gas exposed to a strong external 
radiation field. Here the gas particles are excited to high energy states 
by the absorption of radiant energy. If sufficient collisions could occur, 
this energy would be distributed into an equilibrium distribution of energy 
among the possible states by numerous particle interactions, and the Planck 
distribution and gas absorption coefficient would govern the emission of 
radiation. However, in the rarefied gas, insufficient collisions may occur. 
Emission of energy may then be in a distinctly non-Planckian distribution 
because energy is radiated by spontaneous transitions from these nonequilibrium 
excited states . 

This paper is concerned with the latter case. An approximate grey -gas 
analysis is carried out to delineate roughly the range of variables that 
determines significant departure from LTE. 

ANALYSIS OF NONEQUILIBRIUM RADIATION FROM A GREY GAS 

Consider a volume element dV of grey gas with linear absorption coef- 
ficient a, exposed to radiation from a blackbody source at a temperature T g . 
What is desired is a representation of the spectral distribution of intensity 
emitted by dV if LTE is not assumed. 

If steady state is considered, then the total energy emitted by dV per 
unit time } Q e , must be equal to the energy absorbed by dV per unit time , Q a . 
If dV were in thermodynamic equilibrium, the conservation of radiative energy 
would be expressed by 
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Q a = q e = a(FaTg)dV = 4 a aTg q dV (l) 

where F is the geometric configuration factor between the source and dV, 
and T gq is the equilibrium temperature of the gas . The equilibrium tempera- 
ture of the gas can be related to F and T s by equation (l) as 



If the energy leaving dV per unit time is not in an equilibrium spectral 
distribution, total energy must still be conserved. The total energy leaving 
dV per unit time is then given by 


Q, 


e,ne 


= 4ita dV 



ne 


dA 


(3a) 


Since total energy is conserved, Q e ^ ne must still equal 
tion (l) is 



Q e: ,ne - dV / i A , ne dA = a (Po^) dV 


Q a , which by equa- 
(3b) 


An accurate representation of i A ^ ne can only be obtained by detailed 
analysis of the population of the various energy states in the gas, and the 
various transition processes between these states and their probability of 
occurrence. However, by some judicious assumptions, it is possible to get 
approximate results. 

Consider emission from dV, as an approximation, to be composed of two 
parts. Each part corresponds to equilibrium spectral emission based on a given 
temperature, and these parts will be combined to give the nonequilibrium 
spectral distribution. 

The first part is energy that is absorbed from the incident beam, and 
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is then emitted spontaneously by the gas before being distributed by inter- 
particle collisions- This part is assumed to have a spectral distribution 
similar to that from the blackbody source at T g , which is a Planck distri- 
bution given by i;\ je q^s^* Such an assumption is valid if the spontaneous 
emission causes transition of the particle to the original unexcited state. 
Then energy absorbed at a given wavelength will be spontaneously reemitted 
at the same wavelength unless a collision occurs prior to such emission. 

This criterion of return to the original state is exactly satisfied in the 
two-level atom, which is often used for studies of transition processes, but 
becomes increasingly poor as more and more energy states become possible in 
the gas being considered [5]. 

The remaining fraction of the absorbed energy is assumed to be emitted 
in a spectral distribution ix,eq(^eq) corresponding to the equilibrium tem- 
perature T e( j of dV. This emission is assumed to be from that portion of 
dV that remains in an equilibrium state. Then the intensity at a given 
wavelength is 

i A,ne = A,i A,eq( T s) + ^Ajeq^eq) ( 4 ) 

This synthesis of two equilibrium spectral distributions is analogous to that 
occurring when scattering processes in gray gases are studied. 

Integrating equation (4) over all wavelengths and substituting equa- 
tion (3b) gives 



A'oTo 


ActT 


21 


(5) 


Solving for A' gives, using equation (2) 


6 


A ’ ■ £ - A (%) - (1 - A> (t) (6) 

Equation (2) was derived for LTE. It is valid, however, to use it in the 
non-LTE case, equation (6), because equations (l) and (2) for the relation be- 
tween T g , F, and T eq must also hold for that portion of dV that is assumed 
to be in equilibrium so long as only a gray gas is considered. This is clear 
from examination of equation (l) : If some fraction of Q a is absorbed by the 

equilibrium portion of dV, then the same fraction of Q e in equation (l) must 
be emitted by the equilibrium portion of dV. Equation (2) thus holds for the 
portion of dV assumed to be in thermodynamic equilibrium at temperature 
T g q. Finally, equation (4) can be written 

i A,ne = i A,eq^s^ + Ai A,eq^ T eq^ 

This equation provides the form of i^ ne for the model assumed. To ob- 
tain quantitative results, the value of the parameter A must be estimated. 
This parameter can be interpreted as the probability of a given radiatively 
excited gas particle returning to the equilibrium population based on T eq by 
means of an interparticle collision before undergoing a spontaneous radiative 
emission from the excited state. 

Suppose that there are N 0 excited particles at some given time. If it 
is assumed that some number N R of these excited particles emit spontaneously 
in some time interval t, then assuming a first order rate process, 

dW R = dt (8) 

T R 

where t r is the mean lifetime of a radiatively excited state, and N is 
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the remaining excited population at time t. Similarly, if N c of the 
particles undergo collisions in time t, 

dW c = JL dt (9) 

T c 

where t c is the mean time between collisions. 

Equations (8) and (9) can be integrated to obtain 1% and N c at 
some time t. The probability that the particle deexcites by means of a 
collision is 


A = 



t r + t c 


Substituting this result into equation (7) gives 


L A,ne 


^ 1 + 1 


~fr) ^,eq (T s) + f- i A,eq (T eq ) 

c 


(lo) 


(n) 


Equation (ll) can also be placed in the dimensionless form 



Aj,ne 


T T_ 

^ + (12) 

X eq ? c 


This form has some utility because it illustrates the departure from equilib- 
rium in the emitted spectral distribution of a gas, and is in terms of the 


tabulated blackbody functions 
eters (t r /t g ) and (T s /T eq ). 


X A,eq 


(T)/T 5 


(ref. [6]) and the two param- 
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DISCUSSION OF NONEQUILIBRIUM GREY GAS EFFECTS 

Effect of ratio of transition to collision times . - Examination of equa- 
tion (ll ) shows that the intensity of emitted radiation from a grey gas as a 
function of wavelength is very dependent on the ratio of Tr, the mean life- 
time of an excited state before a radiative transition, to t q, the mean time 
between interparticle collisions. As t r /t c goes to zero, meaning that many 
radiative transitions occur before a collision can take place, the emitted 
spectrum approaches that of the incident radiation reduced in magnitude by 
the factor (T e q/T g )^. As t r /Tq approaches large values, only the last term 
on the right of equation (ll ) is significant, and the emitted spectrum ap- 
proaches that of an equilibrium gas . 

Equation (ll) has the same form as that for the intensity leaving a gas 
with combined absorption and isotropic coherent scattering. For that case 
t r and are replaced by the absorption and scattering coefficients, re- 

spectively [7]. 

In figure 1, the spectral distribution of intensity from a gas volume 
when (T e q/T s ) = 0.1 is plotted for various values of t rA"c> assuming 
T g = 10 000° K. Although the distortion due to the log-log scale must be 
considered, it is obvious that even for T r /t c = 10, the spectrum is quite 
different from the equilibrium Planck spectrum that is present when Tr/t^ -*■ °o. 

The departure from equilibrium is also illustrated in figure 2, again for 
(T e g/T g ) = 0.1. Here, the ratio of i;\,ne i^, e q(^eq) is shown as a func- 

tion of AT e q. For Tr/tq -» 0, which is essentially the case of no collisions 
in the gas being available for the redistribution of energy, the largest devi- 
ation from equilibrium is seen to occur. Emission is in the same spectral 
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distribution as the a&s orbed energy. 

It is interesting that all values of the parameter ^r/^c give an 

intensity ratio of unity at the one AT g q value of 1500 pm °K for this case. 

Thus, if this analysis were completely valid, it would be impossible to 
discern nonequilibrium effects by examining the emitted energy from the gas 
at this AT e q, where T e q would be found in a practical calculation from 
equation (2). Conversely, this is the only AT g q value at which the equilib- 
rium temperature of the gas could be determined by external observation be- 
cause no nonequilibrium effects distort the emitted intensity. 

Also shown on figure 2 is a dotted curve indicating the fraction of 
equilibrium intensity that lies below a given AT e q value. This can be 
used to determine the importance of the departure from equilibrium at various 
AT g q values. For example, the curve for t^/tq = 10 deviates by less than 

10 percent from the equilibrium value for AT e ^ values greater than 1500 pm °K. 

Over 98 percent of the equilibrium intensity lies in this range. Therefore, 
even though the deviation of [i-^ ne A;\ e q(T e q)] from equilibrium reaches 
factors as high as 100 at AT„ n of about 800, the intensity of equilibrium 
radiation at these AT e q values is so small that the spectrum of intensity 
from the gas remains within 10 percent of equilibrium in all important 
spectral regions. 

For t r Ac however, deviations from equilibrium by factors of 10 -3 

occur over AT^ ranges containing significant values of intensity. 

Effect of source temperature . - In the general problems of interest to 
engineers, T e q cannot be specified but rather is defined by equation (2). 

In the discussion presented here, it has been assumed that the interchange 
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factor F Is of such magnitude that it gives the ratios of (T e q/T s ) used in 
the examples. Equations (ll) and (lZ ) could be rewritten to eliminate either 
T g or T e q in terms of the interchange factor, but the present formulation 
seems most useful. In figure 3, the intensity ratio is plotted as a function 
of the source-to-equilibrium gas temperature ratio with Tr/tq as a param- 
eter. This is done for two values of the remaining parameter AT e q. It is 
obvious that the deviation from equilibrium is always greater as 
becomes smaller, regardless of the source-to-gas temperature ratio. The effect 
of (T s /T e q), however, depends upon the specific AT g q value being observed. 

At (T s /T e q) = 1, that is for isothermal conditions, all curves pass 
through i?^ n e/^A, eq(^eq) = 1 so that, as expected, equilibrium is present 
for isothermal conditions regardless of AT or (tr/tcj). 

Parameters governing the departure from equilibrium . - To determine 
readily when the analysis used here predicts significant departure from LTE, 
some additional computations can be made for simple cases. Let i^ eq^eq^ be 
given by Planck's distribution 


i, (T ) = — 
A, eq v eq' , 5 , 


2C, 


2C, 


A [exp(C 2 /AT e q) - 1] A exp(C 2 /AT eg ) 

where the final approximate form, known as Wien's distribution, is within 
1 percent of Planck's distribution for AT less than 3000 qm °K. Substituting 
this final form into equation (13) yields 

*(?0 + ( Tr / Tq ) 


, A,ne 

1 A, eq^ ^eq) 


1 + ( t r/ t c) 


(13) 


where 
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ty(A) - ( qi ] exp 
A s 


AT 


1 - 


eq 1 



(14) 


Use of Wien's distribution should not introduce significant error into equa- 
tion (14) even at AT„ n > 3000 urn °K, because it appears as a ratio. 

For equilibrium to be present, the left-hand side of equation (13) should 
be unity. Examination of equation (13) for this case shows that equilibrium 
is present when either of the conditions 

\|r(A) * 1 (15) 


or 


x T 

rj— » 1 and =2 » i|r(A) 

T c — T C 


(16) 


are met. Condition (15) predicts equilibrium for source and equilibrium gas 

temperatures nearly equal while condition (16 ) shows that for x R sufficiently 

greater than t , other parameters become unimportant. 

The AT, n value at which equilibrium always is present for a given 

(T /T ) is found by substituting equation (14) into equation (15) and solving 
“y. s 

for AT to obtain 


(AT eq ) 


eq 



(17 ) 


Departure from equilibrium in gases under assumed conditions . - For the 
curves of figures 1 to 3, it appears that of IS? or greater would 

assure a close approach to equilibrium energy emission from the grey gas for 
the case of T /t_ of 10. To examine what this means in terms of temperature 
and pressure, a representative case may be examined. 
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Consider a grey gas that attains a temperature of 300° K when exposed 
to a source of radiation at an effective tlackbody temperature of 3000° K. 
First let us see if either of the conditions given by equations (15) and/or 
( 16 ) are satisfied. The value of ^(A) is, using Cg = 14388 pm °K, 



Examination of figure 2 shows that less than 2 percent of the emitted equilib- 
rium energy is at AT e ^ < 1.5X10^, while more than 98 percent lies at 
AT e q < 2X10^ . For these two AT e q values, ty(A) has the values of 0.562 and 
1.91X10“^, respectively. For T p/^Q = 10^, the condition of equation (16) is 
certainly met. What conditions in the gas will give a transition time to 
collision time ratio of this magnitude? 

The value for mean time between collisions in a gas composed of hard 
spheres is approximated from the kinetic theory of gases to be, for a gas 
in equilibrium. 



(19) 


where N is the number of particles per unit volume, cr^ is the mean col- 
lision cross section, and v^ the mean relative velocity of two particles. 
If the number of particles at standard conditions is given by Loschmidt's 
number, Lo = 2.687X10^ atoms per cubic centimeter, the N is given by 
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( 21 ) 


where 

stant, 

as 


v is the mean velocity of the gas particles, k is Boltzmann's con- 
and M the mass of the individual particles. This can be rewritten 



( 22 ) 


where m is the molecular weight of the gas. Substituting equations (22) and 
(20) into (19) gives 


/ VmT 


T C 1.64 OrtP 


-26 

X10 sec 


(23) 


If the gas under consideration were monatomic and of unity atomic weight 
— 1 S 2 

then Cq =» 3X10" cm and m = 1. This provides an upper limit for the 
predicted nonequilibrium effects. For T = 300° K, equation (23) becomes 


r 3.53X10 

T C P 


-11 


sec 


For monatomic gases, is about 10“ seconds [9], so that 


-5. = io 2 


1X10 


-8 


3.53X10 


-11 


or 


P« 0.35 atm 

Thus, for a pure monatomic gas at 300° K exposed to radiation from a blackbody 
source at 3000° K, deviation from equilibrium might begin at as high a pressure 
as l/3 atmosphere, becoming more severe as pressure is decreased. 

For a molecular gas, lifetimes in a given vibrational or rotational energy 
state are on the average much longer than the lifetimes of excited electronic 
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energy states in monatomic gases. Penner [9], for example, quotes molecular 

-1 -3 

energy state lifetimes as being of the order of 10 to 10 seconds. Col- 
lision times for molecules are not significantly different from atoms, how- 
ever, because both m and a ^ increase for molecules, but appear as a ratio 
in equation (23). Therefore, for a molecular gas, deviations from equilibrium 
will not occur until pressures perhaps 10 times smaller than for an atomic 
gas. Such conditions are not often encountered except in studies of radiative 
interaction with interstellar gases. 

CLOSING REMARKS 

The model used in deriving equation (12) has some serious shortcomings 
for quantitative use. Primarily, of course, it is derived for a grey gas, 
a highly unrealistic situation of interest chiefly as an aid to understanding. 
Further, it was assumed that the gas absorption coefficient was unaffected 
by nonequilibrium considerations; in a real gas, this is not true. The 
latter assumption and those outlined in the analysis, however, are not greatly 
in error when the departure from equilibrium is small. Thus, the analysis 
has some utility in predicting the parameter values which signal the departure 
from equilibrium. The equations (15) and (16) thus may give an indication of 
when departure from equilibrium radiative emission might occur. 

LIST OF SYMBOLS 

A*, A emission probabilities defined by equation (4) 

Ci,C 2 constants in the Planck spectral energy distribution 
F geometric interchange factor 

i intensity; energy per unit projected area per unit solid angle per 


unit time 
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i(A) spectral distribution of intensity; intensity per unit wavelength 
k Boltzmann's constant 

Lo Loschmidt's number; 2.687x10 atoms/cm at STP 

M mass of particle 

m molecular weight of particle 

N number of particles per unit volume 

P pressure, atm 

Q energy per unit time 

v mean velocity 

v^ mean relative velocity 

T absolute temperature 

t time 

a linear radiation absorption coefficient 

A wavelength 

a Stefan-Boltzmann constant 

Oq mean collision cross section 

Tp mean time between collision of excited particle with another particle 

mean lifetime of excited state before radiative spontaneous transition 
( A) function defined by equation (14) 

Subscripts : 
a absorbed 

c collisional deexcitation 

e emitted 

eq at equilibrium conditions 

ne at nonequilibrium conditions 
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o at standard conditions 

R radiative deexcitation 

s at source of blackbody radiation 

0 at time t = 0 

A spectrally dependent 
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Fig. 3. - Effect of source temperature on emission from a 
nonequilibrium gas. 


